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ABSTRACT 

This paper provides details of developments 
pertaining to vibration analysis of gyroscopic systems, that 
involves a finite element structural discretization followed 
by the solution of the resulting matrix eigenvalue problem 
by a progressive, accelerated simultaneous iteration 
technique. Thus Coriolis, centrifugal and geometrical 
stiffness matrices are derived for shell and line elements, 
followed by the eigensolution details as well as solution of 
representative problems that demonstrates the efficacy of 
the currently developed numerical procedures and tools. 


NOMENCLATURE 

C c Coriolis acceleration matrix 

C d elastic damping matrix 

K £ elastic stiffness matrix 

K c geometric stiffness matrix 

K' centripetal acceleration matrix 
M elastic mass or inertia matrix 

N shape function 

f initial force vector 

g aeroelastic or structural damping 

«• V-i , imaginary number 

q displacement vector 

q, q velocity and acceleration vectors 

u displacement vector 
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£ damping 

A natural frequencies 

X structural direction cosine matrix 

<f> structural modes of vibration 

n angular velocity matrix, spin rate 

co frequencies 

Superscripts 
T matrix transpose 

-1 matrix inverse 

e elementwise 

INTRODUCTION 

Many practical problems, as turbines, helicopters and 
others are characterized by rotating components arising 
out of functional requirements. Others as spacecraft and 
satellites are often spin stabilized. Their natural frequency 
analysis [1] is of much importance for subsequent 
determination of state of stress as well for stability and 
control analyses. To derive the relevant equations of 
motion one may consider the elastic deformation u ; of a 

point j in a flexible body rotating at a constant angular 
velocity Q about an arbitrary axis, having components 
Q x , Q y , and Q z along the three reference coordinate 

system axes (Fig. 1). Thus, the position vector may be 
derived as 
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whereas the velocity and acceleration vectors can be 
written as 


r)r 

v =— +flxr 

dt 


( 2 ) 


3v 

a = h£2x v 

dt 


d 2 r 

dt 1 

respectively; further 


+ £ixr + 2ftxr + ilx(Qxr) 


(3) 



0 

-O, 


V 

fir = 


0 


<r y 


r a y 


0 

/z. 


(4) 


In Eqn. (3), the first term on the right hand side relates 
to the effects of in-plane stretching on out-of-plane 
deformation; the second term reduces to zero for steady 
spin state while the third and fourth terms relate to Coriolis 
and centripetal acceleration. The preceding equations 
may be conveniently used in formulating strain and kinetic 
energies and the governing equation of motion [2] for the 
entire flexible body may be derived, for the undamped 
case, using Hamilton’s principle as 


Mq + 2MQq + (K e + = -Mfl£2r 


(5) 


or 


Mq + C c q + (K £ + K')q = f c 


( 6 ) 


f c being the centrifugal force vector and r representing 
all Tj, M , K £ , and K' are respectively the inertia, 

elastic and centripetal stiffness matrices, C c being the 
Coriolis matrix. 

The related analysis starts by first computing the 
steady-state deflection due to constant spin £1 by solving 


(K £ +K')q = f c 


(7) 


and thereafter calculating the geometrical stiffness matrix 
K c from resulting in-plane stresses. The equation of 


motion can then be written as below, that also includes 
damping 2 

Mq + (C c +Cjq + (K £ (l + r^)+K , + K G )q = 0 

( 8 ) 

in which C d is viscous damping matrix, g is the structural 
damping parameter, i' being the imaginary number 

V^T. 

To perform relevant dynamic analysis it is essential 
first to derive expressions for nodal centrifugal forces for 
various finite elements spinning about an arbitrary axis. 
Such in-plane forces are used for derivation of k e G 
matrices, where its out-of-plane components are used for 
developing expression for the element centripetal stiffness 
matrix k /e ; examples of some of these matrices are 
derived next. 


NUMBERICAL FORMULATION 
Derivation of Nodal Centrifugal Forces 

Figure 1 shows a typical triangular flat shell element 
undergoing an uniform spin rate Q. R along an arbitrary 
axis. 

Using local coordinate system (LCS) defined by x, y, 
and z axes, the shape function matrix N for in-plane 
motion of the element having components along global 
reference axes X, Y, Z, respectively. 


where 


u = Nu c 


N = RQ 1 


(9) 


( 10 ) 


and in which u is the displacement vector of a typical 
point P within the element; u e is the element nodal 
displacement vector, R is the portion of the shape 
function matrix having elements that are functions of the 
local coordinates x and y , and Q is the portion of the 

shape function matrix with elements expressed in terms 
of the element nodal coordinate values. Further defining 
the internal force vector at any point within the element in 
the LCS as 
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Figure 1 Atypical triangular shell element subjected 
to arbitrary spin rate 


in which l x =cos(X,jc), and so on. Furthermore, 
defining X t , Y t , and Z. as the global coordinates of node 

I of the triangular element, such coordinates for any point 
within the element may simply be obtained as 


X 




X 

Y 

= 

r, 

+x T 

y 

Z 


z,_ 


0 


The in-plane element nodal centrifugal forces may 
next be derived for each of the three components of the 

spin rate Cl R which may then be combined to yield the 
final desired expressions. Thus, typically, corresponding 
the spin component Q . x , the appropriate forces are 
expressed as 


f =l fx fy f J (n) 

, the equivalent concentrated nodal forces of a triangular 
element in the LCS may be expressed as 



P = P&1 


0 

Y 

Z 


(17) 


where p is the mass density. Upon the employment of 

the relationships expressed by Eqns. (15) and (16), Eqn. 
(17) takes the following form: 


which may also be written in the following form: 

f e = | N r fdv 

V 


(13) 


0 


p =p&\ 


m x x + m y y + Y i 
n x x+n y y + Z i 


(18) 


or 


r-jN>* (14) 

V 

where p represents the forces in the global coordinate 
system at any point within an element, having three 
components p x , p Y , and p z in the X , Y , and Z 

directions, respectively, and k is the element direction 
cosine matrix: 



n 

X 

X 



y 

y 

m z 

n z 


(15) 


Expressions for the element nodal forces is then 
derived by appropriate substitution of Eqns. (10), (15), 
and (18) into Eqn. (14), yielding 


f;(n x )=pn 2 x t[Q-'f £ £ R T X dxd y 


(19) 


where t is the element thickness, y k is the y coordinate 

of node k, x t =x k y/y k is the lower bound of the x 
coordinate in the integration process, 
x u =x.~ (x. ~x k )yj y k is the upper bound of the jc 
coordinate, and 


X = 


m x {m x x + m y y+r i )+»,(»,*+»,y+Z/)' 
m y (m x x + m y y + Y,)+n y (n x x+n y y + Z t ) 


( 20 ) 
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(29) 


In Eqn. (19) the vector V p contains the in-plane 
element nodal forces in local x and y directions, being 

expressed as [/„ f y> f Xj f yj 4 f y ] T . 

Associate forces due to spin rates in the Y and Z 
directions may also be similarly expressed as 


fp($ 4 ) = ptfyt [Q-trrw 


( 21 ) 


pQ z t [Q-rn;^ 


(22) 


where 


y = 


X ih x +l y y + X i )+n x (n x x + n y y + Z t ) 
l y (l x x +l y + X i )+n (i n x x + n y + Z ( . ) 


(23) 


$ = 


K (h x + l y y + X,)+m x (m x x + m y y + Z,. )‘ 
l y (l x x + l y + Xj)+m (m x x + m y y + Z i ) 


(24) 


K' e {n z ) = Cl 2 z (l 2 z +m 2 z ) M 

M being the inertia matrix and the complete expression 
for the global centrifugal force matrix is obtained simply as 

K'* = K^(n A .)+K ,c (n i ,)-l-K / ‘ ; (n z ). For quadrilateral 

elements, such data as derived from four constituent 
triangular elements by suitably combining them in an 
orderly fashion. 



The total element in-plane nodal forces are simply 
combined as 


r;=r;(Ox )+»;{«, )+r,'(Oz) (25) 

and may next be transformed into the global coordinate 
system (GCS) as 


P 


= Vf 


T**e 


(26) 


for subsequent stress analysis and generation of the k G 
matrix. 

The out-of plane components of an element 
centrifugal force vector may be derived similarly. Then, 
the element centrifugal force matrices can be readily 
formulated as 


Figure 2 A line element with an arbitrary spin rate 

For a line element (Fig. 2), the coordinate at any point 
within the element is obtained as 


'X' 


X' 


V 

Y 

= 


+x T 

U y 

Z 


X. 


. U z_ 


(30) 


in which u y and u z are deformations at any point on the 
element in the local y and z directions respectively. 

Proceeding as in the case of shell elements and 
assuming nodal lumped masses m , the centrifugal 
forces at a typical i th node expressed in global 
coordinate system. 


K'‘{a x ) = Q. 2 x (n 2 z +m 2 z )M 


(27) 


p‘ x =mCL 2 Y X i + mCl 2 z X i 


(31) 


K' e {Cl y ) = n 2 Y (n 2 z +l 2 z )M 


(28) 


p‘y =mQ. 1 x Y i +m£l 2 7 Y i 


(32) 


p‘ z = mil 2 x Zj +m£l 2 yZ i 


(33) 
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The element centripetal and geometric stiffness 
matrices may then be derived as before. 


(A-/lF)<p = 0 


( 39 ) 


Coriolis Acceleration Matrix 


For any structural element such a matrix can be 
derived as follows 


C e c = 2M e il e 


(34) 


in which F = r*B and can also be expressed in expanded 
form 


Tm 

o' 

-X 

" 0 

* “l> 

Cm 

hr - 

1 ° 

K 


/M 




= 0 


(40) 


in which the element inertia matrix 1VT in calculated by 
employing standard formulation 

Free Vibration of Spinning Bodies 

The associated matrix equation of undamped motion 
of a flexible spinning body, discretized by the finite 
element method, can be derived from Eqn. (8) as 

Kq + Cq + Mq = 0 .... 


where F is a pure Hermitian matrix, the roots X = i* co 
are real and occur in pairs X l , -X l , ■■■, X^, -X N , 

associated eigenvectors occurring in complex conjugate 
pairs. Eigenvalues of the original Eqn. (38) may simply be 

obtained as Xj i , while noting that eigenvectors are the 
same for both the cases. 

Solution of Eqn. (39) has been achieved earlier by 
various solution techniques (3,4,5]. A variant of the multi- 
vector progressive iteration technique [6,7,8] is presented 
next. 


in which K = K £ +K , + K C is the stiffness matrix, 
C = C c the Coriolis force matrix, and M is the inertia 
matrix. For small vibration the matrices K and M are 
real, symmetric and positive definite, whereas C c is 

skew symmetric, being also a function of the spin fit. 
Eqn. (35) may be resolved into a set of first-order ordinary 
differential equation as below resulting in a linear 
eigenvalue problem. 


Ay + By = 0 


(36) 


in which 



"M 

O' 


0 

-M' 

A = 



, B = 



0 

K 


M 

C 



The matrices A and B are of order (2 xjV), N 

being order of K , M, and C matrices, noting that B is 
skew-symmetric. Solution of Eqn. (36) may be assumed 
as y = ye 0 ” yielding 


(A + oM) (p-O ^ 

where the eigenvalues (ti are pure imaginary, the 
eigenvectors occurring in complex conjugate pairs. Eqn. 
(38) may be further rearranged [3] as 


Eiqensolution by a Progressive multivector Iteration 

Equations (38) and (39) are used for the computation 
of the first NR roots and vectors. Associated solution 
steps are given next: 

Step 1. Perform Cholesky factorization of the real 
symmetric and positive definite matrix A of order 
2Nx2N, taking into account the sparsity of the 
constituent matrices, initially only, 


A = LDL r 


(41) 


L is an unit lower triangular matrix, D being diagonal 
matrix. 


Step 2. Form a set of randomly generated 
2x NRT (NRT > NR) real, trial vectors each of order 
2 N , NRT being number of trial vectors. 


where 


1 /V 2 A 

y x 


NRT * 


NRT 


y] 


(42) 


J m—J 


9 


x+i* j y, i" =^^A 


(43) 


Step 3. Initially form 
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(44) 



NRT 


i,l 


Step 4. Using only unconverged roots, solve 


AY /+I =BX, 


(45) 


or in individual matrix form 


KY® = MX) u) + CX)° 


(46) 


yielding Y$ 


MY = -MX) 0 


(47) 


that is = -X)° . 

Similarly solve 


AX i+I =BY,. +1 


(48) 


which can be elaborated as below 


KX^MY^+CY/'} 


(49) 


Step 6. Setting ‘p^x, 2 p=y, 3 p= 2 x, 4 p= 2 y,---, 
orthonormalize each y p i+1 vector with respect to the A 
matrix, only for unconverged ( NRl-NR ) vectors, 
yielding p, as follows: 


2 xnri r 

'Pt+i = J Pm ~ £ [('P J 


Pm 


P- 


/=1 


U A'p (+1 ]'p (+I 


l=2xNR\+l 


(52) 


; P I+ i = %i/[('P i+ iF A'Pm 


1/2 


j = 2xNRl + 1, 2xNRT 


Then set back ‘x='p, ‘y= 2 p, 2 x= 3 p, 2 y= 4 p, and 
so on. 

Step 7. Go to step 4 to continue iteration till solution 
convergence is achieved, in which 

V NRl+2* NRT* 

A i+1 _ L X i+1 X i+ 1 ■” X i+lJ (53) 

Step 8. The Rayleigh-Ritz method, involving only the 
unconverged roots, is employed next to solve the further 
reduced eigenvalue problem 


yielding X)J.j and also 


MX 


(«) _ 
i+l — 


-MY 


(0 

i+l 


(50) ' n which 


A; + iQ, + i — F i+1 Q i+1 A 


(54) 


from which X)“} = -Y^j 

In above the superscripts u and l refer the upper 
and lower half of a vector, respectively. 


A /+1 = 


F , = <£> , 

i+i /+i i+i 


(55) 


with 


Step 5. Using unconverged vectors only, perform 
vector norm check by computing an estimate of the 
magnitude of a typical vector as 



and check root convergence if 
{ J A i + l - J A i )/ j A m < EPSN , a prescribed convergence 
parameter. Assuming that first NR\ roots have 
converged, if A7?l = NR , proceed to step 8. 



|w«i+i 


<P,+i 


NR 1+2* 

9 ,+, 



(56) 


A being an approximation of A , the diagonal matrix 
containing the eigenvalues. 

Step 9. A root convergence test is then performed on 
each unconverged root under consideration, 




/ 


j JL 


2+1 


< EPS 


(57) 


in which EPS is the root convergence factor. 
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Step 10. Recompute the eigenvectors, 


1. Spinning Cantilever Beam 


® i+1 ®/+iQ;+i (58) 

and relabel <I> (+1 as 4> i+1 for possible use in step 4 

/V A 

involving X ;+1 components of <J> 1+1 . 

Step 11. Goto step 4 and perform steps 4, 6, 8-10 
only if the number of converged roots NRl < NR . 
Otherwise go to step 12. 

Step 12. End of analysis. 

The analysis procedure described herein enables 
effective computation of the first few roots and vectors 
which are often the requirement for the important 
subsequent stability and dynamic response analyses. For 
further references, steps 4-7, that involves checking root 
norm convergence, will be referred to as phase 1 
operation. Subsequent phase 2 operations, involving 
Rayleigh-Ritz root convergence calculations consist of 
step 4, 6, and 8-11. Thus, progressive saving in solution 
time is obtained in the following key areas, (a) back 
substitution operation in step 4, (b) vector norm check 
process in step 5, (c) mass orthonomalization calculation 
in step 6, (d) Rayleigh Triple matrix multiplication in step 
8, (e) Rayleigh vector computation in step 9. It may be 
noted that all major computations in steps 1-7 are 
performed in real number thereby saving considerable 
computational efforts. Solution accuracy factors for 
EPSN , EPS , and EPSR are set to 0.05, 0.00001 and 
0.0001, respectively. 


NUMBERICAL EXAMPLES 

Three example problems are presented next that 
testifies to the applicability and efficacy of the currently 
developed solution algorithms for computation of roots 
and modes of rotating bodies. The first example relates 
to a spinning cantilever beam subjected to an arbitrary 
spin rate. Results are presented for both the undamped 
beam as well as solutions in the presence of viscous and 
structural damping. The second and third example 
problems pertain to a rotating cantilever plate, involving 
rather large degrees of freedom. Solution results and 
CPU times are compared for the usual and the 
accelerated version of the associated code. 


Figure 3 shows a cantilever beam spinning about an 
arbitrary axis. Associated data parameters, in consistent 
British units, are as follows: 


Young’s modulus, E = 30x1 0 6 

Cross-sectional area, A =1.0 

Moment of inertia, 

about y axis = 1/12 

about z axis = 1/24 

Mass density, p = 1/6 

Length, L = 60.0 

Scalar nodal viscous damping = 0.6283 

(in x , y , z translations) 

Structural damping coefficient = 0.01 

Spin rate along Y axis(£2 K ) = 0.1 Hz 


Y 

A?’ 

2 3 4 5 6 7 8 9 10 1 1 

z / s 

Figure 3 Spinning cantilever beam 

Figure 4 depicts the first few mode shapes. 


Y 



Figure 4 First few mode shapes for a cantilever beam 

(£ly = 0 ) 


Table 1 lists the natural frequencies of the cantilever 
beam with various forms of damping. 
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Table 1 Natural frequencies of a spinning cantilever 
beam 


0) 

"O 

o 

Natural frequencies, A (rad/sec) 

*« 

II 

o 

P 

ll 

0.6283 (rad/sec) 


undamped undamped 

Viscous damping 

1 

2.6631 ±2.7626 

- 0.3799 ± / 2.7365 

2 

3.7660 ±3.7932 

- 0.3800 ±/ 3.7742 

3 

16.5039 ±16.8409 

- 0.3872 ±h 6.6360 

4 

23.3387 ± 23.7476 

- 0.3870 ± /23.7441 

5 

45.7488 ±47.0130 

- 0.3937 ±/47.01 10 

6 

64.6946 ± 66.3958 

- 0.3933 ± / 66.3944 

7 

88.7152 ±92.0851 

- 0.4012 ±/ 92.0840 

8 

125.4547 ± 130.0557 

- 0.4005 ±/*1 30.0550 

9 

145.0221 ±152.3365 

- 0.4105 ±h 52.3357 

10 

205.0799 ±215.0823 

-0.4093 ±/21 5.081 7 


Natural frequencies, A (rad/sec) 

<D 

"O 

o 

£2 y = 0.6283 (rad/sec) 


Structural damping 

Viscous and 
Structural damping 

1 

0.0138 + / 2.7626 

- 0.3659 + / 2.7365 

-0. 0138-/2.7626 

- 0.3938 - Z2.7365 

o 

0.0189 + / 3.7932 

- 0.3608 + / 3.7742 

c. 

-0. 0189-/3.7932 

- 0.3990 - / 3.7742 

3 

0.0842 + /*1 6.8411 
-0. 0842 -/*1 6.8411 

- 0.3030 + / 1 6.6361 

- 0.4714 -M 6.6361 

4 

0.1187 + / 23.7479 
-0. 1187-/23.7479 

— 0.2682 + / 23.7444 
-0.5057-/23.7444 

5 

0.2351 + / 47.01 35 
-0. 2351 -/ 47.0 135 

-0.1585 + / 47.01 15 
-0.6287-/47.0115 

6 

0.3320 + / 66.3966 
- 0. 3320 - / 66.3966 

-0.6130 + / 66.3952 
- 0.7252 - Z66.3952 

7 

0.4604 + /92.0862 
- 0. 4604 - / 92.0862 

- 0.5922 + / 92.0850 
-0.8616-/92.0850 

8 

0.6502 + /*1 30.0574 
-0. 6502 -h 30.0574 

- 0.2497 + /‘l 30.0565 

- 1.0508 -h 30.0565 

9 

0.7617 + M 52.3384 
-0. 7617 -M 52.3384 

- 0.3512 + /'l 52.3375 
-1.1721 -h 52.3375 

10 

1.0754 + / 21 5.0850 
-1.0754-/215.0850 

-0.3799 + /21 5.0843 
-1.4847-/215.0843 


It may be noted that for the undamped spinning cases 
as well as in the presence of structural damping the roots 
occur in pairs as A l ,-A l ,A 2 ,~A 2f --- 7 A, N ,-Z N , the 

latter being also complex in nature. For the viscous 
damping case the roots are in complex conjugate form 

whereas when the viscous and structural 
damping are combined, the real part tends to be 
somewhat different corresponding to identical imaginary 
parts, while also noting that the roots of the original set of 

Eqn. (30) are defined as 0) = A./C . 

Note: lump mass option for non-spinning cases, 
consistent mass option for spinning cases. 


2. Spinning Plate 1 

A square plate (Fig. 5), cantilevered at the base is 
discretized by a 20x20 mesh, and analyses are conducted 
for spin rates along four sets of local axes (£2 X , 

£l z , and Q R ). The spin rate along an arbitrary axis, £l R 
(at 45° of each axis) = 100 rad/sec, with an input of Q. x = 

Cl Y = Cl z = 57.735 rad/sec; also 100 rad/sec for 
individual axis spin rates along x , y , and z local axes. . 

Associated data parameters, in consistent British units, 
are as follows: 


Young’s modulus, E 
Side length, L 
Thickness, t 
Poisson’s ratio, jU 
Mass density, p 


10x10 6 

10.0 

0.1 

0.3 

0.259x1 O' 3 


Y 



X 


Figure 5 Spinning cantilever plate 
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Results of these analyses are summarized in Table 2. 

The first 12 frequencies are given in the table; first few 

mode shapes are also shown in Fig. 6. 3, Spinning Plate 2 



(a) Mode 1 (b) Mode 2 



(c) Mode 3 (d) Mode 4 



(e) Mode 5 (f) Mode 6 

Figure 6 First 6 mode shapes for a cantilever 
plate (Q, = 0) 


Table 2 Natural frequencies (A. in rad/sec) of a 
spinning cantilever plate (20 x 20) 
(consistent mass option) 


In an effort to compare the efficacy of the accelerated 
progressive, process over the original simultaneous 
iterative technique, another plate (Fig. 7) that results in a 
rather larger order problem is analyzed next. Associate 
details of the plate, expressed in British units, are as 
follows: 


Side length, L 
FE mesh 

Spin rate, (f2 z ) 


= 30.0 

= 30 x 30 

= 100 rad/sec 


Y 



X 


Figure 7 Cantilever plate (30x30) spinning along Z- 
axis 


CD 

■Q 

Non 

Cl z =100 

£l R = 100 

£l x =100 

Sly = 100 

O 

spinning 






(ft = 0) 





1 

216.80 

245.20 

162.02 

126.19 

51.73 

2 

529.15 

549.78 

515.09 

497.29 

496.24 

3 

1325.91 

1350.84 

1331.52 

1333.44 

1309.70 

4 

1689.95 

1703.87 

1690.48 

1681.95 

1685.67 

5 

1922.04 

1943.08 

1928.36 

1928.88 

1913.21 

6 

3352.53 

3367.97 

3358.46 

3357.42 

3350.10 

7 

3810.79 

3835.60 

3823.53 

3830.02 

3804.62 

8 

3970.11 

3981.71 

3974.04 

3968.21 

3972.17 

9 

4401.49 

4422.87 

4412.28 

4415.73 

4398.49 

10 

5732.01 

5744.32 

5737.70 

5735.50 

5732.96 

11 

5993.17 

6012.35 

6003.40 

6006.02 

5992.21 

12 

7375.74 

7397.47 

7388.79 

7393.43 

7374.15 


Other details are as in the previous plate example 
problem. The results in a 5400 degrees of freedom 
problem and the relevant eigensolutions are effected 
using the current progressive simultaneous iterative 
solution method as well as the original simultaneous 
iterative solution; this affords a quantitative comparison of 
the two analysis procedures. Each eigensolution is 
derived to compute the first 100 modes and roots using 
110 randomly generated trial vectors. For both solutions, 
8 steps were needed to achieve convergence in phase 1 
and another 1 0 steps were needed in phase 2 to achieve 
a final solution convergence. Table 3 provides a list of 
first few frequencies, for both spinning and non-spinning 
cases. 
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Table 3 First 20 natural frequencies (A. in rad/sec) of a 
30x30 spinning plate (consistent mass option) 


Mode 

D 

II 

o 

Q z = 100 

1 

24.1007 

109.2439 

2 

58.9329 

131.7314 

3 

147.6217 

210.4748 

4 

188.3865 

310.7209 

5 

214.3126 

347.9099 

6 

374.5417 

414.7329 

7 

424.7105 

519.2821 

8 

443.5959 

573.6855 

9 

491.2739 

647.7871 

10 

641.6125 

734.3848 

11 

670.5685 

816.0489 

12 

824.2067 

902.8594 

13 

858.9417 

961.7461 

14 

890.9186 

997.2900 

15 

962.1829 

1145.9511 

16 

1033.7367 

1180.3106 

17 

1089.4643 

1276.6296 

18 

1358.6855 

1387.4848 

19 

1377.8232 

1440.7086 

20 

1381.0980 

1476.7978 


Analyses were conducted for two combinations of 
progressive features; Table 4 presents relevant details. 
These data indicates that the progressive method is more 
efficient than the usual simultaneous iteration technique. 


Table 4 Solution details for a 30x30 spinning plate (T: 
true, F: false for progressive feature) 


Computation Phase 

PSI 1 

SI 2 

Phase 1 

T 

F 

Phase 2 

T 

F 

Iteration Steps in Phase 1 

7 

8 

Iteration Steps in Phase 2 

11 

11 

CPU time (sec) 

78.17 

141.79 


CONCLUDING REMARKS 

The paper presents the novel derivations of 
constituent Coriolis, centrifugal and geometrical stiffness 
matrices for various finite elements undergoing rotation 
along arbitrary axes. This is followed by the matrix 
formulation of the associated equation of motion, which 
proves to be convenient for subsequent natural frequency 
analysis. An accelerated simultaneous iteration solution 
method is also presented in some details. Finally, a 
number of example problems are presented for a beam 
and plates rotating along a set of arbitrary axes, that 
testifies to the efficacy and broad applicability of the 
currently developed techniques. A comparison of CPU 
time is also provided for the current accelerated solution 
against the original simultaneous iterative method, noting 
the marked improvement. These techniques were 
incorporated in a multidisciplinary finite element analysis 
code [2], capable of solving a large array of practical 
problems. 
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